Spin fluctuations, susceptibility and the dipole oscillation of a nearly ferromagnetic 

Fermi gas 
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We discuss the spin fluctuations and the role played by the magnetic susceptibility in an atomic 
Fermi gas interacting with positive scattering length. Both thermal and quantum fluctuations are 
considered. Using a sum rule approach and recent ab initio Monte Carlo results for the magnetic 
susceptibility of uniform matter we provide explicit predictions for the frequency of the spin dipole 
oscillation of a gas trapped by a harmonic potential and discuss the deviations from the behaviour 
of an ideal gas when the system approaches the ferromagnetic transition. The role of the Landau's 
parameters in the characterization of the magnetic properties is also discussed. 
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The recent experiment on itinerant ferromagnetism 
reported in Ref. 1] has stimulated a novel debate on 
the magnetic properties of interacting Fermi gases [2- 
|8[. In this experiment it vifas shown that, by increasing 
the value of the positive scattering length a, the system 
exhibits characteristic features of the Stoner instability 
where the gas enters a new phase of ferromagnetic type. 
The experiment was performed on a metastable state of 
a Fermi gas. Actually, due to the attractive nature of 
the interaction, there exists a 2-body bound state for 
positive values of the scattering length and the ground 
state of the system is a Bose gas of molecules (see, e.g., 
[9[ and references therein). Indeed, three-body collisions, 
which convert atoms into molecules, provide one of the 
signatures of the transition to a magnetic phase, since 
they become more relevant by increasing a in the nor- 
mal state, but are less efficient in the magnetic phase 
due to the formation of spatial domains of non interact- 
ing atoms. They give rise to a maximum in the atomic 
losses. Other signatures of the phase transition were a 
maximum in the cloud size and a minimum in the kinetic 
energy [l| . This experiment did not however measure di- 
rectly any magnetic property of the system. In particular 
no evidence for magnetic domains was found above the 
transition point. Recent Monte Carlo simulations [ld.lll| 
have shown that the ferromagnetic transition of an inter- 
acting Fermi gas takes place at the value kpa ^ 0.8 — 0.9, 
where kp = {Sir'^nY'^ is the Fermi momentum and n the 
density of the sample. This value is significantly smaller 
than the experimental value found at the lowest temper- 
atures. This discrepancy, together with the absence of 
evidence for magnetic domains, suggests that the physi- 
cal behaviour of the system is far from being understood. 
In particular the occurrence of 3-body losses [j, |7| and 
the competition between the tendency of the system to 
become ferromagnetic and to form pairs [8J are likely at 
the origin of new physical effects. 

One of the most relevant thermodynamic quantities 
characterising the transition to the ferromagnetic phase is 
the magnetic susceptibility x ~ {Vd'^A/d{N^ — A'^j.)^)"^, 



where A is the free energy, V the volume of the sample 
and A't and A^i are the number of fermions in the two 



H 



spin-states. In general the interaction will increase the 
susceptibility of the gas (see, e.g., [12|) which, in a second 
order phase transition, should diverge at the transition. 
The actual order of the transition is sensitive to temper- 
ature, disorder and polarisation (see, e.g., [a Isl Il3l Il4| 
and reference therein) 15]. 

At finite temperature the spin susceptibility can be 
related to the measurement of the thermal spin fluctua- 
tions: 



AiN,-N^ = k,T^ 



N 



(1) 



where N — N^ + N^ and n is the density of the sample. 
The measurement of the spin fiuctuations, and through 
(HI) of X, seems to be now accessible, as it has been shown 
by the very recent shot noise measurements of the den- 
sity fluctuations in an ideal Fermi gas Ig, |l7[ . In such 
experiments the fluctuations are measured in a small sub- 
volume of the whole atomic cloud, where the system can 
be considered uniform [18[. 

At low temperature the applicability of ([T]) becomes 
less and less efficient because of the suppression caused 
by the thermal factor and the emergence of quantum fluc- 
tuations which exhibit a weaker dependence on N, but 
become soon important as T — > 0. The quantum fluctua- 
tions of the number of particles have been investigated in 
[l9|. The result for the spin fiuctuations, to the leading 
order in N , can be written as 
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where the parameter a is fixed by the low q behavior of 
the static structure factor, according to S°'{q — >■ 0) = 
aq/kp and C is a cut-off constant determined by the 
short range behaviour of S°'{q). In the case of the ideal 
Fermi gas one finds a = 3/2'^/3, C ~ 10.45. Thus the 
comparison with Eq. ([T]) using the T — ideal Fermi gas 
susceptibility value xo = 3n/2ep, with ep the Fermi en- 
ergy, shows that quantum fluctuations become dominant 



for temperatures T < Tf/N^'^{().^2 + 0.12 InA''), where 
Tp = tp/kB is the Fermi temperature. Thus, taking 
for example, T ~ Q.2Tp quantum fluctuations become 
important already for N < lO-^. 

For interacting Fermi gases the value of x ^^'^ ^ 
can be calculated employing Landau theory of Fermi 
liquids. Moreover the coefficient a is bounded by 
the spin susceptibility through the sum rule inequality 
a < \/epx{T = 0)/2n. Thus, an important question is 
whether and how also a - and consequently the T = 
quantum fluctuations - diverges at the transition point 
identified by the divergence in x- 

In Landau theory of Fermi liquids the susceptibility is 
given by 
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where m* = m{l + F^ /3) > m is the effective mass, and 
Ff, F° are the l-th angular momentum symmetric (den- 
sity) and antisymmetric (spin) Landau parameters, re- 
spectively. The susceptibility and thus the thermal fiuc- 
tuations ([T]) present a divergence for Fq = — 1 [20|. 

The parameter a can also be expressed in terms of 
the Landau parameters. Indeed, in general, S°'{q) = 
— 7r~^ J /mx(g, w) where x(?)'^) is the linear response 
function determined by the complete quasi-particle scat- 
tering amplitude [21|. It is well known that the linear 
resp onse function presents a pole for a; = at F^ = — 1 
|2l| . By taking into account only the / = Landau pa- 
rameters, a reads 
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and, thus it also diverges for Fq = — 1, but just logarith- 
mically, a ex ln(l/(l -|- Fq)). This means that quantum 
fluctuations increase more slowly than thermal fluctua- 
tions as one approaches the critical point of the ferro- 
magnetic transition. 

When the interaction is small the Landau parameters 
can be determined using perturbation theory |22J and, to 
second order in kpa, they read j23i] 



F^ = -ikpa) + -^i2 + ln2){kFa)^, (5) 

F^^^{7\n2-l){kFa)\ (6) 

l{kpa)-^{l-ln2){kpa)\ (7) 

F^ = ~^{2 + ln2){kpa)\ (8) 
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so that the inverse susceptibilty, to the same order, takes 
the form 
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(9) 
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FIG. 1: Inverse spin susceptibility for an homogeneous Fermi 
gas as a function of the interaction parameter kpa. Points: 
Monte-Carlo data from [ly]. Dashed (green) line: fit to the 
Monte-Carlo data. The dot-dashed (blue) and the continuous 
(red) lines are the first and second order expansion, respec- 
tively (see Eq. ([9])). 



To the same order one can also calculate the inverse com- 
pressibility of the gas: k'^/k = (1 -|- F^)m/m* , where 
kq = Z/2nep the ideal Fermi gas compressibility. For 
higher values of kpa the susceptibility can be calculated 
numerically using quantum Monte-Carlo techniques 'lC|. 
The results for x~^ ^■^ ^ function of the dimensionless pa- 
rameter kpa are shown in Fig. [T] where the Monte Carlo 
predictions are compared with the small kpa expansion 
^. The figure shows that, a ferromagnetic instability 
appears around {kpa)c — 0.83. The first order expan- 
sion of x~^i extrapolated to the critical point, gives the 
well known Stoner result {kpa)c = 7r/2 ~ 1.57, while the 
second order expansion would give {kpa)c — 1.05. 

Aside from the fluctuations, another way to investi- 
gate the magnetic properties of Fermi gases is the study 
of their dynamic behaviour as a response to a spin per- 
turbation. In this letter we analyse the zero-temperature 
spin dipole oscillation of a two component Fermi gas in 
the presence of an harmonic trap ^(r) = l/2m{uj'^x'^ + 
ojyu'^ + ojIz^). The dipole mode for a weakly interacting 
Fermi gas has been already studied in [24,] and experi- 
mental results were reported in [25|. In this letter we 
estimated the frequency of this oscillation using a sum 
rule approach [26| based on the calculation of the ratio 



^SD 



mi 



(10) 



where m^ — ^^ |(0|Z3|n)p(i?„ — £'0)'^ arc the moment 
of the strength distribution function relative to the spin 
dipole operator D — Ylii'^i'X ~ ^4)- The ratio (jTOl) ac- 
tually provides an upper bound to the frequency of the 
lowest energy state excited by the operator D. Differ- 
ently from uniform matter, where a negative value of the 



spin parameter Fq gives rise to Landau damping, in the 
presence of harmonic trapping one expects a discretized 
cohective mode with frequency smaher than the free os- 
cillator value. 

The energy weighted sum rule entering the numerator 
of pUj) is easily calculated in terms of a double commu- 
tator: 



mi 



[D,[H,Dm^N 



2m 



(11) 



where we have used the fact that only the kinetic energy 
operator gives a contribution to the commutator [H, D] . 
The inverse energy weighted sum rule is instead related 
to the spin dipole polarizability and can be calculated, 
at zero temperature, by minimizing the total energy of 
the system in the presence of an external static coupling 
of the form —XD. In the local density approximation the 
energy can be written as: 

E= [ dr[eK(r), n;(r)) - Xz{n^{r) - n^{r))] (12) 
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FIG. 2: Main panel: Spin dipole frequency as a function of 
the interaction parameter iPpa (see text). The dashed line is 
the first order expansion Eq. {2)). Inset: the value of the 
interaction parameter kpa in the trap as a function of fc^a. 
kp is the Fermi momentum of the interacting cloud calculated 
at the centre of the trap. 



where t{n^^^n\^ is the energy density of uniform matter. 
By expanding e(n-f , n^) up to quadratic terms in n^ — nj,, 
minimization of E yields the result n-|- — nj, = \z chi(n) 
for the polarization induced by the external field where 
X~"'^(n) = d^t/din^ — TT'iY' is the zero temperature in- 
verse magnetic susceptibility of uniform matter calcu- 
lated at the local value of the density. The calcu- 
lation of the induced dipole spin moment then yields 
m_i = 1/2 J drz^xii^) a-nd hence the result 



N 



^SD 



m J drz'^xi'^) 
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for the dipole frequency. Eq. (J13p shows explicitly that 
an increase of the mag netic susceptibility will result in 



a decrease of ujsd [27[. The density profile, n(r), en- 



tering the integral P^ should be calculated in the local 
density approximation using the equilibrium condition 
/i(n) -I- Vho = Mo where ^{n) = de{n)/dn is the chemical 
potential and /Xq is fixed by the normalization condition. 
One can easily check that, using the ideal gas expres- 
sion for the magnetic susceptibility, Eq. P^ yields the 
expected result ujsd = ^z where Wz is the trapping oscil- 
lator frequency in the z-th direction. Inclusion of interac- 
tion effects modifies the value of the spin dipole frequency 
both through the change of the magnetic susceptibility x 
and the change of the density profile fixed by the equa- 
tion of state n{n). At the first order in the scattering 
lenght we get 
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where aho is the geometrical average of harmonic oscil- 
lator lengths. Result p^ coincides with the prediction 



made in [24| employing a scaling ansatz. Such a coinci- 
dence is not surprising since one can easily show that, in 
the limit of the ideal Fermi gas, the polarization induced 
by the static dipole field coincides with the scaling form 
n^ — n^ — {X/muj'^)dzn. According to perturbation the- 
ory result ^^ then corresponds to an exact result for 
the spin dipole frequency calculated to first order in a. 

In Fig. [5] we report the value of the spin dipole fre- 
quency (fT5)) calculated using the Monte-Carlo results for 
^J,{n) and x{^) oi uniform matter [10| and using the local 
density approximation to account for the inhomogene- 
ity caused by the external potential. The results are re- 
ported as a function of kpa, where kp = (2ANY/^ /auo is 
the Fermi momentum of an ideal Fermi gas calculated in 
the centre of the trap, aho is the the geometrical average 
of the harmonic oscillator lengths. In the inset we report 
the value of the interaction parameter kpo. as a function 
of k'pa. The instability, which appears initially in the 
centre of the trap at around fc^a ~ 1 corresponds to the 
(Monte-Carlo) bulk value kpa ~ 0.83 calculated in the 
centre of the trap. The dipole frequency is seen to be 
significantly affected by the interaction even for values of 
/c^a well below the ferromagnetic transition. This result 
is promising in view of the possibility of checking experi- 
mentally our understanding of the magnetic properties of 
the Fermi gas in the regime where the expansion in pow- 
ers of k'pa should be universal (i.e. independent of the 
actual form of the two-body potential) and where the life 
time of the sample is large enough [ij . At the transition 
point the dipole frequency (|10|) does not vanish since the 
divergence of the magnetic susceptibility concerns only 
the behaviour of the gas in the centre of the trap where 
the integral (fT5|) is suppressed by the factor z^. 



It is worth noting that the result p3|) provides an up- 
per bound to the lowest spin dipole frequency. An even 
lower bound can be obtained by noticing that, accord- 
ing to Landau theory of Fermi liquids, the /-sum rule in 
the spin channel is not exhausted by the low-lying spin 
excitations. Indeed the low-energy quasi-particle excita- 
tions give a contribution to the /-sum rule proportional 
to the combination m/m*(l -|- l/SF^^"''). Thus, while in 
the density channel the sum rule is not affected by the 
interaction, the spin dipole /-sum rule is instead reduced 
by the interaction since Ff > Ff , the remaining part of 
the /-sum rule being provided by multi-pair excitations 
located at higher energy [28| . The renormalized value of 
the mi sum rule up to the second order in kpa can be 
explicitly determined using the expressions ^ and ^ 
for the Landau parameters. We find: 
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-(l + 81n2)(fcF(r)a)2 



(15) 
The explicit calculation of the integral ([T5]) yields a re- 
duction of the spin dipole frequency of less than 10 per 
cent at fc^a = 0.5 with respect to the prediction of (fT3|). 

Let us conclude by pointing out that an open issue 
concerns the effects of collisions on the damping of the 
spin dipole mode and in particular its dependence on 
temperature and on the value of the scattering length 
expecially close to the transition [29[. 
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